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Abstract.  The  interaction  of  a  gravity  wave  with  a  solar  tide  is  analyzed  using  ray  theory 
in  order  to  assess  whether  the  temporal  oscillation  of  the  tide  has  any  significant  effects  on 
the  interaction.  We  consider  two  types  of  solution:  a  “full  ray  solution,”  in  which  tidal 
accelerations  are  included  in  the  gravity  wave  ray-tracing  equations,  and  a  second 
“Lindzen  solution,”  in  which  they  are  neglected;  the  latter  is  so  named  because  it  yields 
similar  results  to  the  parameterization  of  Lindzen  [1981].  Initially,  we  consider  an  idealized 
tide  of  constant  velocity  amplitude  in  a  steady  isothermal  atmosphere,  as  this  allows 
analytical  solutions  to  be  derived.  A  numerical  ray-tracing  code  is  employed  to  determine 
ray  solutions  within  more  complex  tidal  backgrounds.  Full  ray  solutions  often  differ 
markedly  in  amplitude,  wavenumber,  and  trajectory  from  Lindzen  solutions  for  the  same 
wave,  highlighting  the  importance  of  tidal  accelerations.  Tidal  accelerations  have  a 
stabilizing  influence  on  gravity  wave  amplitudes  by  refracting  waves  to  larger  intrinsic 
phase  speeds,  thus  reducing  both  the  occurrence  and  intensity  of  tidally  modulated  gravity 
wave  breaking.  Even  in  the  absence  of  dissipation,  time-varying  refraction  gives  rise  to 
time-oscillating  gravity  wave  action  densities,  HA /lit,  which  in  turn  lead  to  wave 
momentum-flux  densities  which  are  anticorrelated  with  tidal  winds  in  agreement  with 
mesospheric  observations.  Although  dA/dt  averages  to  zero  over  one  cycle  of  a  constant- 
amplitude  tide,  its  phase  locking  to  the  tide  implies  transient  flow  accelerations  which 
temporarily  amplify  tidal  wind  oscillations  and  might  permanently  change  tidal  amplitudes 
in  a  nonlinear  interactive  model.  Numerical  simulations  through  more  complex  tidal 
backgrounds  show  instances  where  permanent  dA/dt  values  are  induced  in  a  gravity  wave 
through  its  interaction  with  the  tide.  Our  results  point  to  weaknesses  in  a  Lindzen 
parameterization  of  gravity  wave-tidal  interactions.  A  simple  extension  of  the 
parameterization  is  suggested  which  incorporates  some  of  these  time-varying  influences. 


Introduction 

Dissipation  of  gravity  wave  amplitudes  in  the  atmosphere 
results  in  an  Eliassen  and  Palm  [1960]  (EP)  flux  divergence 
which  either  accelerates  or  decelerates  the  local  flow.  In  the 
middle  atmosphere  this  wave  driving  of  the  flow  tends  to  occur 
almost  continuously,  which  maintains  global  flow  patterns  that 
differ  substantially  from  a  purely  radiative  circulation.  The 
dissipated  wave  energy  is  transferred  into  turbulence  which 
diffuses  chemical  species,  which  in  turn  influences  the  global 
and  seasonal  distributions  of  trace  gases  in  the  middle  atmo¬ 
sphere. 

Consequently,  mean-flow  accelerations  and  turbulent  diffu¬ 
sion  produced  by  dissipating  gravity  waves  need  to  be  incor¬ 
porated  accurately  into  global  models  of  the  middle  atmo¬ 
sphere.  Lindzen  [1981]  parameterized  the  physics  of  these 
processes  by  assuming  a  single  linear  saturated  wave.  His  pa¬ 
rameterization  proved  amenable  to  numerical  implementa¬ 
tion,  and  those  models  of  the  middle  atmosphere  which  incor¬ 
porated  it  often  exhibited  major  improvements  in  their 
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simulations  [see,  e.g.,  Garcia  and  Solomon,  1985;  Palmer  et  al., 
1986],  These  successes,  along  with  its  simplicity  and  general 
physical  basis,  quickly  made  it  the  standard  way  of  incorporat¬ 
ing  gravity  wave  processes  into  large-scale  models.  Despite 
this,  various  shortcomings  in  the  parameterization  have  been 
noted.  These  include  its  omission  of  localized  turbulence  ef¬ 
fects  [e.g.,  Chao  and  Schoeberl,  1984;  Coy  and  Fritts,  1988; 
Walterscheid  and  Schubert,  1990],  supersaturation  [e.g., 
Lindzen ,  1988;  McIntyre ,  1989],  transience  [e.g.,  Fritts  and 
Dunkerton,  1984;  Walterscheid  and  Schubert,  1990;  Laprise, 
1993],  and  various  nonlinear  and  nonhydrostatic  effects  [e.g., 
Fritts,  1985;  Walterscheid  and  Schubert,  1990;  Keller,  1994], 
Recently,  parameterizations  of  a  dissipating  spectrum  of 
gravity  waves  have  been  developed  [Fritts  and  Lu,  1993;  Med¬ 
vedev  and  Klaassen,  1995].  While  these  spectra  are  more  typical 
of  the  average  wave  energy  distributions  encountered  in  the 
middle  atmosphere,  individual  wave  events  still  often  dominate 
the  EP  flux  divergence  in  the  stratosphere  [e.g.,  Bacmeister  et 
al.,  1994]  and  (less  frequently)  in  the  mesosphere  [e.g., 
Yamamoto  et  al.,  1987;  Fritts  et  al.,  1988].  Thus  variants  of  the 
Lindzen  scheme  are  still  used  commonly  within  middle  atmo¬ 
sphere  models  to  generate  realistic  diffusivities  and  circulation 
patterns  [see,  e.g.,  Garcia  and  Boville,  1994;  Jackson  and  Gray, 
1994;  Huang  and  Smith,  1995;  Kinnersley,  1996], 

In  the  mesosphere  the  wave-breaking  situation  is  further 
complicated  by  the  presence  of  large-amplitude  solar  tides. 
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Dissipation  of  tidal  amplitudes  can  also  produce  significant 
mean-flow  accelerations  in  the  equatorial  mesosphere  and 
lower  thermosphere  [e.g..  Miyahara  et  al.,  1993;  Hamilton , 
1995],  Furthermore,  observations  have  revealed  that  large- 
amplitude  tidal  wind  oscillations  modulate  gravity  wave  ampli¬ 
tudes  and  EP  flux  divergences  [Ganguly,  1980;  Frills  and  Vin¬ 
cent,  1987;  Wang  and  Fritts,  1991;  Thayaparan  et  al.,  1995],  It  is 
aspects  of  this  latter  process  that  wc  focus  upon  here. 

The  way  in  which  gravity  waves  might  interact  with  tides  was 
explored  initially  by  Spizzichino  [1970]  and  Wallerscheid  [1981], 
Frills  and  Vincent  [1987]  applied  a  simple  gravity  wave  satura¬ 
tion  parameterization  to  model  the  interactions  that  they  ob¬ 
served  in  radar  winds,  which  revealed  diurnally  varying  EP  flux 
divergences  in  the  presence  of  a  diurnal  tide.  Bjarnason  et  al. 
[1987]  also  applied  the  Lindzen  parameterization  in  a  photo¬ 
chemical  model  to  simulate  tidally  modulated  diffusividcs  due 
to  gravity  wave  breaking.  These  diurnally  varying  diffusividcs  in 
turn  produced  a  diurnal  variation  in  mesospheric  ozone  abun¬ 
dances  which  agreed  well  with  observations.  Later  studies  us¬ 
ing  the  Lindzen  parameterization  illustrated  that  the  mere 
presence  of  tides  significantly  altered  even  the  daily  averaged 
mean-flow  accelerations  and  diffusivitics  produced  by  gravity 
wave  breaking  [Hunt,  1986,  1990;  McLandress  and  Ward,  1994], 

Fritts  and  Vincent  [  1 987]  also  modeled  the  effect  that  tidally 
modulated  gravity  wave  dissipation  might  have  upon  the  tide 
itself.  Reduced  tidal  amplitudes  and  an  advancement  of  tidal 
phase  in  time  were  predicted.  A  similar  gravity-wave-induced 
damping  of  mesospheric  Rossby  waves  was  postulated  by 
Lindzen  [1984]  and  Schoeberl  and  Strobel  [1984]  and  later  sim¬ 
ulated  in  a  model  by  Miyahara  [1985]  using  the  Lindzen  pa¬ 
rameterization  (see  also  Miyahara  etui.  [1986]).  Tidal  modelers 
have  used  Lindzen-likc  parameterizations  to  examine  the  ef¬ 
fects  of  gravity  wave  dissipation  on  global  tidal  structures  [e.g., 
Forbes  et  al.,  1991;  McLandress  and  Ward,  1994].  Lu  and  Fritts 
[1993]  used  a  spectral  parameterization  to  assess  the  effects 
(see  also  Roble  and  Ridley  [1994]):  They  too  found  advances  in 
tidal  phase  but  noted  that  increases  and  decreases  in  tidal 
amplitudes  were  possible  within  their  model.  Progress  in  the 
modeling  of  these  interactions  has  been  reviewed  by  Miyahara 
and  Forbes  [1994],  Fritts  [1995],  and  Walterscheid  [1995], 

Use  of  the  Lindzen  parameterization  in  the  presence  of  a 
tide  requires  that  the  tidal  oscillation,  as  “seen”  by  a  gravity 
wave,  can  be  validly  approximated  as  vertically  varying  but 
time-independent,  so  that  conventional  nontransient  wave  for¬ 
mulae  can  be  used  to  describe  the  propagation  and  amplitude 
of  any  gravity  wave  through  the  tide  [e.g.,  Lindzen,  1981].  This 
implies  that  the  tidal  phase  and  group  speeds  must  be  much 
slower  than  the  group  velocity  of  the  gravity  wave.  Here  we 
assess  the  accuracy  of  these  assumptions  by  computing  gravity 
wave  trajectories  in  the  presence  of  tides  using  ray  techniques. 
We  use  ray-tracing  equations  appropriate  for  non  hydrostatic 
gravity  waves,  as  given  by  Marks  and  Eckermann  [1995],  which 
we  extend  here  to  incorporate  the  refractive  effects  of  tidal 
winds  which  oscillate  both  spatially  and  temporally.  Through¬ 
out  this  study  we  simplify  the  assumed  mean  and  tidal  back¬ 
ground  wherever  possible  in  order  to  focus  on  the  essence  of 
the  gravity  wave  response  within  these  tidal  wind  fields. 

In  drafting  this  manuscript  the  results  of  a  similarly  moti¬ 
vated  ray-tracing  study  by  Zhong  et  al.  [1995]  came  to  our 
attention.  Their  numerical  simulations  showed  that  time- 
varying  tidal  oscillations  produced  important  refractive  effects 
on  gravity  wave  trajectories,  as  we  have  also  found.  Given  this, 
wc  have  changed  the  focus  of  this  paper  somewhat  to  concen¬ 


trate  less  on  ray  trajectories  (many  examples  of  which  are  given 
by  Zhang  et  at.  [1995])  and  more  on  the  effects  on  gravity  wave 
amplitudes  and  gravity  wave-tidal  interactions,  which  Zhong  et 
al.  [1995]  did  not  simulate.  Additionally,  we  have  attempted  to 
present  a  simple  physical  picture  of  the  processes  which  oper¬ 
ate  here  by  deriving  analytical  solutions  for  a  gravity  wave  in 
the  presence  of  an  idealized  tidal  wave  of  constant  amplitude, 
drawing  on  studies  of  related  phenomena  in  the  ocean.  In  this 
way,  comparisons  among  these  results  and  those  from  simpler 
nontransient  wave  formulae,  such  as  the  Lindzen  parameter¬ 
ization,  are  more  easily  made. 

Simplified  Ray-Tracing  Model 

We  employ  ray  techniques  to  model  gravity  wave  (GW) 
propagation  through  a  prescribed  tidal  oscillation.  General 
simulations  of  this  kind  can  be  performed  accurately  with  the 
Gravity-wave  Regional  or  Global  Ray  Tracer  (GROGRAT), 
the  first  version  of  which  (V1.0)  was  described  by  Marks  and 
Eckermann  [1995].  GROGRAT  is  a  general  numerical  algo¬ 
rithm  which  simulates  gravity  wave  propagation  through  an 
arbitrary  atmosphere  using  ray-tracing  methods.  We  are  pres¬ 
ently  engaged  in  a  major  upgrade  of  GROGRAT  to  “next- 
generation”  status  (V2.lt  and  beyond).  One  purpose  of  this 
upgrade  is  to  include  time  variations  of  the  background  atmo¬ 
sphere  realistically  into  the  model,  to  complement  the  three- 
dimensional  spatial  variations  which  are  already  catered  for. 
The  present  study  provides  a  simplified  platform  for  develop¬ 
ing  and  testing  time-varying  ray  equations,  while  simulta¬ 
neously  providing  useful  insights  into  a  specific  phenomenon 
of  considerable  current  interest  in  the  middle  and  upper  at¬ 
mosphere. 

While  general  numerical  ray-tracing  models  like  GRO¬ 
GRAT  and  others  [e.g.,  Zhong  et  al.,  1995]  give  accurate  solu¬ 
tions,  they  arrive  at  these  solutions  by  numerically  accumulat¬ 
ing  influences  from  a  range  of  background  parameters  and  all 
of  their  spatiotempora!  derivatives.  Thus  it  is  not  always  easy  to 
identify  and  separate  the  major  dynamical  interactions  from  a 
host  of  other  minor  contributions  to  these  ray  solutions,  yet 
this  is  important  in  potentially  refining  the  simple  paramctcr- 
izations  of  these  processes  used  in  larger-scale  models.  Thus  to 
aid  theoretical  interpretation,  we  shall  use  a  “stripped-down” 
version  of  the  GROGRAT  ray-tracing  code,  in  which  we  retain 
only  those  processes  which  are  essential  to  a  first-order  de¬ 
scription  of  the  GW-tidal  interaction. 

Ray  Equations 

The  ray  equations  used  here  are  as  given  in  Appendix  A  of 
Marks  and  Eckermann  [1995].  For  example,  refraction  of  the 
vertical  wavenumber  m  is  given  by  the  ray-tracing  equation 


■  (i) 

which  replicates  equation  (A3f)  of  Marks  and  Eckermann 
[1995]  and  uses  identical  notation:  k  =  ( k ,  l,  m)  is  the 
wavenumber  vector,  «+  =  w  —  k  •  U  is  the  intrinsic  frequency, 
and  U  is  the  background  flow  velocity,  currently  given  an  as¬ 
sumed  horizontal  form  ( U ,  V,  0).  Additional  background 
variables  arc  the  Brunt-Viiisala  frequency  A,  the  Coriolis  pa¬ 
rameter  /,  and  a  =  I/27/p,  where  H  is  the  density  scale 
height.  The  time  derivative  in  the  ray-following  frame  is  d/dt, 
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while  subscripts  denote  partial  derivatives  in  the  Earth-based 
frame  (i.e„  Uz  —  dU/dz),  and  A  =  k2  +  l2  +  m 2  +  a2. 

When  local  time  variations  of  these  background  quantities 
(e.g.,  U,  -  aUiiH)  are  ignored,  the  ground-based  wave  fre¬ 
quency  o)  remains  constant  along  the  ray  ( du/dl  =  0)  [see, 
e.g.,  Lindzen,  1981;  Marks  and  Eckermann,  1995].  However,  a 
time-varying  background  atmosphere  due  to  the  presence  of 
tides  results  in  an  additional  ray-tracing  equation  [e.g.,  Jones, 
1969;  Whitham,  1974] 

da)  1 

Jf^kU'  +  IV,  +  +  l2)  -  <a2),(*>+2  -  f2)]- 

(2) 

Thus  dw/dt  #  0  in  general,  and  u>  is  no  longer  constant  along 
the  ray  path. 

Use  of  ray  methods  to  specify  gravity  wave  paths  requires 
that  the  background  atmosphere  varies  slowly  compared  to  the 
wave  oscillation,  so  that  the  (WKB)  assumption  of  a  single 
slowly  modulated  wave  packet  remains  valid.  Thus  following 
Marks  and  Eckermann  [1995]  (see  also  Broutman  [1984]),  we 
first  require  that 


inequalities  required  of  both  our  exact  calculation  of  (3)  and 
our  approximate  calculation  of  (4). 

Our  provision  for  a  time-varying  background  atmosphere 
means  that  we  must  also  satisfy  slowly  varying  requirements  in 
the  time  domain,  which  yields  the  following  temporal  analogue 
of  (3)  [e.g.,  Broutman,  1984]: 

1  I  I 


Note  that  there  is  no  temporal  analogue  of  the  Einaudi  and 
Hines  [1970]  WKB  derivation  which  led  to  (4).  Equation  (5) 
can  also  be  evaluated  exactly  using  dm/dz  values,  as  shown  in 
Appendix  A. 

Simplified  Mean  and  Tidal  Background 

In  order  to  study  the  essence  of  the  GW-tidal  interaction  we 
consider  an  idealized  zonal  wind  field  of  the  form 

U(x,  z,  t)  =  £/t[de(z)  sin  (Kx  +  Mz  +  nClt  +  <p), 

~  ^tide(y,  z,  t ),  (6) 


Since  we  evaluate  dm/dz  directly  in  this  study,  as  described  in 
Appendix  A,  then  the  along-ray  calculation  used  by  Marks  and 
Eckermann  [1995]  to  approximate  (3)  is  not  needed  here. 
Zhong  et  al.  [1995]  used  the  more  formal  WKB  parameter  of 
Einaudi  and  Hines  [1970], 

S2  =  <2i  +  Qz  «  1.  (4) 

where  Q,  —  -0.15m'  4(dm/dz)2  and  Q2  —  0.5 m~2  d2m/dz2. 
While  Q  |  resembles  S2  in  (3),  Q-,  is  an  additional  higher-order 
term.  Zhong  et  al.  [1995]  found  that  Q2  was  nonnegligible  in 
their  simulations,  thus  implying  that  the  ray  method  of  solution 
may  be  questionable  in  many  circumstances.  We  estimated  Q2 
in  several  similar  simulations  by  subtracting  dm/dz  at  adjacent 
heights  along  the  ray  path  to  estimate  an  approximate  <)2m/ 
dz2  value  and  also  found  that  nonnegligible  Q  ,  values  some¬ 
times  arose. 

There  are  two  points  worth  noting  about  using  along-ray 
calculations  of  (4)  to  gauge  the  applicability  of  the  ray  method. 
First,  the  extra  terms  in  (4)  account  for  possible  partial  reflec¬ 
tion  of  wave  energy  due  to  sudden  changes  in  the  “refractive 
index”  for  the  waves  [see,  e.g.,  Blumen,  1985],  Significant  par¬ 
tial  reflection  of  wave  energy  does  not  necessarily  invalidate 
the  use  of  ray  methods  to  continue  tracking  the  transmitted 
wave  energy  (although  energy  correction  terms  should  obvi¬ 
ously  be  included).  A  second  point  is  that  evaluating  Q2  by 
differencing  adjacent  points  along  the  ray  path  does  not  give 
exact  vertical  partial  derivatives,  since  such  a  calculation  im¬ 
plicitly  includes  horizontal  and  time  variations  in  m  as  well 
(assuming  a  generally  oblique,  time-varying  ray  path).  As  men¬ 
tioned  in  Appendix  A,  accurate  calculations  of  spatiotemporal 
partial  derivatives  of  wave  parameters  at  various  points  along 
the  ray  path  generally  require  some  sort  of  “ray  tube”  tech¬ 
nique,  in  which  the  parameters  of  a  number  of  closely  spaced 
rays  are  interpolated  to  form  a  “tube,”  from  which  these  de¬ 
rivatives  can  be  inferred  numerically  [e.g.,  Marks  and  Ecker¬ 
mann,  1995],  So  as  not  to  complicate  things  unnecessarily  in 
this  Study,  we  use  (3)  as  our  guide  (as  in  section  5  of  Einaudi 
and  Hines  [1970])  but  ensure  that  our  calculations  obey  the 


where  K  —  n/(a  cos  0)  is  the  tidal  horizontal  wavenumber  (a 
is  the  Earth’s  radius,  8  is  latitude,  and  n  is  a  positive  integer), 
M  is  the  tidal  vertical  wavenumber,  and  Q  =  277  (24  hours)-1. 
Thus  the  flow  is  produced  by  a  purely  zonal  tidal  wind  oscil¬ 
lation  of  peak  amplitude  CTrDE(z). 

We  take  the  background  density  p0  at  z  —  60  km  to  be  3  X 
10-4  kg  m-3  and  adopt  a  steady  isothermal  background  tem¬ 
perature  T  =  240  K,  following  McLandress  and  Ward  [1994], 
Variations  in  T  (and  hence  N  and  a)  can  be  large  at  the 
equator,  where  they  make  important  contributions  to  gravity 
wave-tidal  interactions  [e.g.,  Miyahara  and  Forbes,  1994].  How¬ 
ever,  comparisons  of  isothermal  and  anisothermal  extratropi- 
cal  simulations  by  Zhong  et  al.  [1995]  usually  revealed  small 
differences  in  ray  trajectories,  so  we  omit  such  effects  here  for 
simplicity  and  clarity  of  subsequent  interpretation.  One  should 
note,  however,  that  these  and  other  parameter  dependences 
(e.g.,  mean  wind  variations)  can  be  accurately  gauged  using 
GROGRAT  in  the  presence  of  more  complicated  background 
atmospheres. 

For  the  simplified  flow  (6)  then,  since  K  is  wavenumber-/? 
(i.e,,  of  planetary  scale),  Ux  =  dU/dx  is  small,  so  k  remains 
nearly  constant,  while  l  is  refracted  only  slightly  due  to  the 
latitudinal  variation  of  /  [e.g.,  Dunkerton,  1984;  Marks  and 
Eckermann,  1995],  Thus  xh  —  (k2  +  l2)U2  remains  constant 
with  height  to  a  very  good  approximation  when  tracing  any 
small-scale  gravity  wave  through  these  tidal  winds. 

Consequently,  to  simplify  things  further  in  this  study,  we 
have  set  K  =  0  in  all  the  simulations  to  be  reported  here,  so 
that  the  flow  becomes 

U(z,  t)  =  UTlDE(z)  sin  (Mz  +  nklt  +  <p), 

=  £2tide(^i  t)-  (7) 

We  have  conducted  simulations  with  and  without  the  Kx  term 
in  (6)  and  have  found  only  small  differences  in  the  solutions 
when  the  waves  do  not  propagate  over  long  horizontal  dis¬ 
tances.  The  neglect  of  horizontal  tidal  variations  is  also  ad¬ 
dressed  more  formally  in  deriving  analytical  solutions  for  a 
constant  tidal  amplitude  in  the  following  section. 

An  advantage  of  using  (7)  is  that  it  permits  direct  compar¬ 
isons  of  our  results  with  those  of  earlier  studies,  which  used 
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zonal  flow  fields  identical  to  (7)  but  did  not  incorporate  the 
effects  of  time-varying  tidal  winds  on  the  gravity  wave  equa¬ 
tions  fe.g.,  Walterscheid,  198!;  Lu  and  Fritts,  1993;  McLandress 
and  Ward ,  1994].  Thus  any  differences  which  arise  can  now  be 
directly  ascribed  to  time-varying  tidal  effects  acting  on  the 
gravity  wave.  Another  advantage  is  that  the  dependence  of  the 
flow  onz  and  l  alone  allows  accurate  determinations  of  gravity 
wave  amplitudes  to  be  made  along  the  ray,  as  described  in  the 
following  section. 


Gravity  Wave  Amplitudes 

By  far  the  greatest  difficulty  involved  in  including  back¬ 
ground  time  variations  into  a  ray-tracing  approach  involves  the 
accurate  calculation  of  wave  amplitudes  along  the  ray  path.  A 
solution  to  the  wave-action  continuity  relation 


c)A 

~di 


+  V  ■  (cs/4)  = 


2A 

T 


(8) 


is  required  [Bretherton  and  Garrett,  1969],  where  A  is  the  action 
density,  cB  =  (c,/x,  cf/y,  c,/z)  is  the  ground-based  group  velocity 
vector,  and  r  1  is  a  damping  term.  Within  the  simplified  back¬ 
ground  environments  that  we  have  just  introduced,  horizontal 
gradients  of  the  horizontal  action  flux  terms  in  (8)  disappear, 
so  that  (8)  can  be  simplified  and  reexpressed  as 


dA 

~di 


(9) 


where  cyz  =  dzldt  =  dw/dtn.  Again,  dcf)z/dz  in  (9)  can  be 
evaluated  generally  using  a  ray  tube  calculation.  However, 
within  the  restricted  backgrounds  considered  here,  and  on 
manipulating  the  equations,  this  term  can  be  evaluated  accu¬ 
rately  along  single  ray  paths  using  a  ray  equation  for  dm/dz 
and  given  a  sufficiently  accurate  initial  condition.  The  method 
is  outlined  in  Appendix  A  and  is  used  in  this  model. 

The  damping  rate  r_l  in  GROGRAT  parameterizes  the 
effects  of  background  turbulence  and  radiative  damping  on 
wave  amplitudes  [sec  Marks  and  Eckennann,  1995,  Appendix 
B],  but  in  the  interests  of  simplicity  we  omit  such  influences 
here.  However,  we  must  include  wave-breaking  effects  in  some 
way.  To  compare  as  closely  as  possible  with  the  Lindzcn 
scheme,  we  replace  the  nonhydrostatic  saturation  scheme  used 
in  GROGRAT  with  the  simpler  static-instability  saturation 
limit  used  by  Lindzen  [1981],  which  yields  «sat  =  1  using  the 
nomenclature  of  Marks  and  Eckennann  [1995].  The  saturated 
wave-action  density  then  becomes  [ Marks  and  Eckermann, 
1995] 


A 


Sill 


\ _Po(z)w1  r  p_  Al2+(0  +  2/i  f2 
4  (k*  +  P)  [ 1  +  co+2  +  A2- tu  1  2  {  1 


(10) 

Thus  whenever  A  >  Hsat,  action  density  is  dissipated  by  in¬ 
serting  a  damping  factor  t  =  rsaI  into  (9)  which  returns  the 
wave  to  marginal  stability  (A  =  ASM). 

Types  of  Solution 

In  what  follows  we  shall  consider  two  main  solutions.  The 
first  is  the  solution  produced  using  the  ray  methods  just  de¬ 
scribed  and  will  be  referred  to  as  the  full  ray  solution. 

The  second  is  calculated  by  keeping  w  constant  (dw/dt  =  0), 
whereupon  m  simply  follows  from  the  dispersion  relation 


m  ==  K'h  a/GTp--  (n) 

Group  velocities  are  evaluated  from  these  wave  parameters 
and  are  used  to  define  the  ray’s  trajectory  in  this  solution. 
Wave  amplitudes  are  calculated  by  solving  the  wave-action 
equation  in  the  absence  of  time  variations,  which  yields  the 
solution 


Acv  =  const.  (12) 

This  second  ray  solution  gives  results  which  are  similar  in  many 
respects  to  the  parameterization  formulae  of  Lindzen  [1981], 
as  Schoeberl  [  1 985 1  demonstrated.  Hence  we  shall  refer  to 
them  as  “Lindzen  solutions,”  although  one  should  bear  in  mind 
that  the  correspondence  is  not  exact.  For  example,  the  nonhy¬ 
drostatic  dispersion  relation  (11)  represents  a  slight  extension 
of  the  hydrostatic  Lindzen  theory. 

In  a  steady  environment  with  purely  vertical  background 
wind  variations  the  full  ray  solution  and  the  Lindzen  solution 
give  identical  results.  A  major  goal  of  this  study  is  to  assess 
vvhether  the  differences  between  the  two  solutions  that  arise  on 
adding  lime  variations  to  the  background  are  significant.  We 
shall  focus  most  closely  on  any  differences  in  wave  amplitudes 
and  the  predicted  amounts  and  locations  of  wave  breaking. 


Analytical  Ray  Solutions  for  Constant 
Tidal  Amplitude 

We  consider  first  the  simplest  case  of  CT1DB(z)  equal  to  a 
constant,  f)T1DE,  in  (6).  In  this  case  the  ray  and  amplitude 
equations  for  the  gravity  wave  can  be  solved  analytically.  So¬ 
lutions  lo  a  somewhat  similar  oceanic  problem  were  derived  by 
Thorpe  [1989]  and  tested  using  tank  experiments. 

The  geometry  of  this  problem  is  depicted  in  Figure  1.  It 
shows  a  larger-scale  gravity  wave  (in  our  ease,  a  solar  tide)  of 
constant  velocity  amplitude  propagating  to  the  right,  with 
phase  moving  down  (Cz  <  0)  and  energy  moving  upward 
(C,/z  >  0).  We  assume  a  transverse  wave  with  nonzero  K,  so 
we  must  also  incorporate  a  vertical  tidal  velocity  perturbation 
W’  (see  Figure  1).  A  smaller-scale  gravity  wave  impinges  upon 
the  larger-scale  wave  oscillation  from  the  bottom  left.  Assum¬ 
ing  both  waves  are  aligned  east-west,  then  from  (8)  the  action 
density  A  of  the  small-scale  wave  must  obey  the  continuity 
relation 


r)A  d  d 

+  +  iiZ  =  (13) 


in  the  absence  of  dissipation.  AH  three  terms  in  (13)  are  clearly 
nonzero  on  inspecting  Figure  1.  However,  on  rotating  the  axes 
through  an  angle  [3  =  arctan  (  —  KIM),  such  that  the  new  axes 
x'  and  z'  are  aligned  parallel  and  orthogonal,  respectively,  to 
the  phase  fronts  of  the  tidal  oscillation  (see  Figure  1),  (hen  (13) 
transforms  and  simplifies  to 


dA 

<11 


+  dz~  ~  0, 


(14) 


where  cgz,  =  c!/2  cos  /3  -  cgx  sin  J3  is  the  gravity-wave  group 
velocity  along  z'.  The  derivation  of  (14)  follows  from  the 
invariance  of  A  under  Galilean  transformations  fe.g.,  Grirn- 
shaw,  1984]  and  since  d/dx'  (cgx,A)  =  0,  as  can  be  noted  on 
inspecting  Figure  1. 

Finally,  we  transform  to  a  third  nonaccelerated  frame  with 
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axes  (x",  z” ).  These  axes  have  the  same  tilted  orientation  as 
the  (x',  z' )  axes  in  Figure  1  but  move  downward  along  the  z’ 
axis  at  the  tidal  phase  speed,  so  that  x"  =  x'  and  z"  =  z‘  — 
Cz,l ,  where  Cz,  =  —  |C|  (see  Figure  1).  Then  (14)  transforms 
and  simplifies  to 

(cm„A )  =  0,  (15) 

since  tidal  time  variations  disappear  in  this  frame  (i.e., 
dA/tit  =  0).  Equation  (15)  has  the  solution 

c,)t„A  =  [(«■„  cos  fi  -  cm  sin  fi)  -  C.,]  A  =  const, 

(16) 

where  Cj(  =  nCl/HC,  and3£  =  ( Kz  +  M2)1'2.  Note  that  cgz  = 
W'  +  c nz  and  cux  =  U'  +  cgx,  where  (cgx,  c~z)  is  the  intrinsic 
group  velocity.  However,  since  W'  cos  fi  —  U'  sin  fi  —  0  (see 
Figure  1),  then  (16)  can  also  be  expressed  as 

((cj,  cos  /3  -  tv  sin  |3)  -  C.,}  A  =  const,  (17) 

This  same  transformation  sequence  also  simplifies  the  ray¬ 
tracing  equations  of  the  small-scale  wave,  such  that  in  the 
frame  (x",  z")  the  observed  frequency  to"  and  wavenumber 
component  alongjt”,  denoted  k",  remain  constant.  The  follow¬ 
ing  Doppler  relations  hold  among  the  frequencies  of  the  small- 
scale  gravity  wave  measured  from  these  different  reference 
frames: 

to  =  to'  =  to"  +  C,,m' , 

=  to"  +  Cz,(m  cos  fi  -  k  sin  fi).  (18) 

Since  K2  «.  M2  for  solar  tides,  then  fi  =  arctan  (  —  KIM)  is 
very  small  and  the  rotational  transformation  (x,  z )  —*  (x',  z') 
produces  a  very  small  effect  in  this  problem.  This  implies  that 
the  tidal  gradients  along  the  x  axis  are  small  and  have  little 
influence  on  the  ray  path  and  amplitude  of  the  small-scale 
gravity  wave. 

Thus  we  omit  these  effects  by  setting  fi  =  0,  just  as  we  did 
earlier  when  setting  K  =  0  in  going  from  (6)  to  (7)  and 
assuming  W  =  0.  This  simplifies  the  transformation  sequence 
to  a  single  transformation  from  (x,  z)  to  (x",  z"),  such  that 
z"  —  z  —  Cj ,  x"  —  x,  and  the  axes  (x" ,  z")  move  downward 
with  the  vertical  phase  velocity  of  the  tide  Cz  -  nVUM .  The 
solutions  now  simplify  to  those  derived  and  analyzed  by  Bronf¬ 
man  and  Young  [1986]  in  their  analysis  of  an  oceanic  gravity 
wave  propagating  through  an  inertial  oscillation.  Again,  the 
wave  frequency  as  measured  from  the  frame  (x",  z")  remains 
constant,  so  that  from  (18),  fi  =  0  yields 

u>"  =  (»  —  C;  m  =  const.  ( 19) 

Thus  the  ground-based  frequency  a>  is  periodically  refracted  as 
the  tide  oscillates  temporally,  according  to  (19),  which  is  the 
solution  of  (2)  for  this  problem.  The  remaining  wave  parame¬ 
ters  vary  so  as  to  maintain  the  dispersion  relation 

io+ 2  —  [w"  +  Cz  m  —  k  U(z,  r)]2 

(. k 2  +  l2)N2  +  ( m 2  +  a2) / 2 
=  - . — - £ - —  •  (20) 

Using  (16)  with  fi  =  0,  the  wave-action  solution  now  simpli¬ 
fies  to  [Broutman  and  Young,  1986] 


Figure  1.  The  contour  plot  depicts  the  velocity  oscillations  of 
a  large-scale  gravity  wave  with  a  peak  velocity  amplitude  of 
1 1  m  s  ' 1  (contour  labels  are  in  meters  per  second  and  negative 
values  are  shaded).  The  wavenumber  vector  3C,  velocity  oscil¬ 
lation  vector  ({/',  0,  W'),  group  velocity  C,/r  and  phase  ve¬ 
locity  C  for  this  wave  are  also  depicted.  A  smallcr-scale  gravity 
wave  of  wavenumber  k  impinges  on  the  big  wave  from  the 
bottom  left,  with  phase  and  group  velocities  of  c  and  crj,  re¬ 
spectively,  as  shown.  The  rotated  axes  (x',  z')  are  used  to 
derive  amplitude  solutions  for  the  small  wave  as  it  propagates 
through  the  big  wave. 

A(Cp  -  Cr)  =  const,  (21) 

where  cgz  —  cgz.  This  solution  determines  the  amplitude  of 
the  small-scale  gravity  wave  along  its  path  through  the  tide  in 
the  absence  of  saturation.  These  solutions  provide  a  conve¬ 
nient  initial  check  of  our  numerical  ray  method  of  evaluating/! 
within  more  general  environments,  as  described  in  Appendix  A. 

Note  that  (21)  is  singular  whenever  cl/z  -  C ,  (i.e.,  A  -»  oo). 
This  "phase-group  condition”  produces  ray  caustics  and  has 
been  interpreted  in  oceanic  studies  as  a  reflection  point  in  the 
(x",  z")  frame  which  produces  “separatrix  crossings”  and  per¬ 
manent  changes  in  m  and  co+  [Broutman  and  Young,  1986; 
Broutman,  1986;  Broutman  and  Grimshaw,  1988;  Bruhwiler  and 
Kaper,  1995],  Phase-group  caustics  do  not  occur  here  because 
both  the  tide  and  the  gravity  wave  always  propagate  their 
energy  upward,  so  c,)Z  and  C\  are  antisigned  and  never  equal. 

Results  for  Constant  Tidal  Amplitudes 

Diurnal  Tide 

We  consider  first  a  diurnal  tide  ( n  =  1)  of  peak  amplitude 
0Tide  =  20ms'1  and  Az  =  2tt/M  =  40  km.  These  tidal 
parameters  are  fairly  typical  of  the  midlatitude  winter  meso¬ 
sphere  [e.g.,  Vincent  et  al.,  1988;  Manson  et  al.,  1989],  so.  we 
choose  y  -  45°N  initially.  Figure  2a  plots  contours  of  the 
height-time  zonal  velocity  oscillation  of  the  tide. 

The  remaining  curves  throughout  Figure  2  show  the  height 
variations  of  a  range  of  gravity-wave-related  variables  which 
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Figure  2.  Variations  with  height  along  the  ray  path  of  various  wave  parameters  of  a  gravity  wave  propagating 
through  tidal  zonal  velocities  (contour  labels  in  Figure  2a  are  in  meters  per  second  and  solid  (dotted)  contours 
depict  positive  (negative)  tidal  winds).  The  results  of  two  ray  solutions  are  depicted:  the  Lindzen  solution 
(dashed  curve)  and  the  full  ray  solution  (solid  curve).  At  z  =  60  km,  u 77  +  x/7  =  7  m2  s~2  for  both  solutions, 
(a)  Tidal  velocities  are  contoured;  curves  depict  ray  loci  through  these  diurnal  tidal  winds,  (b)  spatial 
trajectory,  (c)  intrinsic  frequency,  (d)  vertical  wavelength,  (e)  ground-based  zonal  phase  speed,  (f)  intrinsic 
zonal  phase  speed,  (g)  unsaturated  action  density,  (h)  saturation-scaled  action  density,  (i)  Eliassen-Palm  flux 
divergence,  and  (j)  mean-flow  acceleration. 


resulted  from  tracing  a  gravity  wave  through  these  tidal  winds. 
The  solid  curves  show  the  full  ray  solution,  while  the  dashed 
curves  show  the  Lindzen  solution  for  the  same  wave.  In  each 
instance  wc  launched  the  wave  zonally  (/  =  0)  at  z  =  60  km 
and  t  —  0  with  an  initial  horizontal  wavenumber  k  —  2ir 
(40  km)-1,  an  initial  ground-based  horizontal  phase  speed 
cx  —  to/k  of  25  m  s_l,  and  an  initial  horizontal-velocity  vari¬ 
ance  u77+  x/7  =  7  m2  s-2.  The  wave  is  launched  at  a  trough 
of  the  tidal  velocity  oscillation,  at  which  dmldz  -  0  initially 
(see  Appendix  A). 

The  resulting  lime-height  ray  loci  are  superimposed  upon 
the  tidal  velocity  contours  in  Figure  2a,  while  the  spatial  tra¬ 
jectories  of  each  solution  are  shown  in  Figure  2b.  Note  that  the 


group  trajectories  of  the  two  solutions  differ  considerably. 
These  differences  reflect  the  important  effects  of  tidally  in¬ 
duced  background  accelerations  dU'riV)K/dl  on  the  ray  path  of 
the  gravity  wave.  In  other  words,  the  Lindzen  solutions  give 
inaccurate  ray  paths  through  the  tide. 

The  remaining  panels  in  Figure  2  show  the  variation  with 
height  along  the  ray  path  of  various  gravity  wave  parameters. 
We  note  that  the  refractive  excursions  in  o>+  and  A,  =  2 Trim 
(Figures  2c  and  2d,  respectively)  are  smaller  for  the  full  ray 
solution  than  for  the  Lindzen  solution.  The  reason  for  this  is 
illustrated  by  the  cx  =  10/k  variations  in  Figure  2e.  While  the 
ground-based  zonal  phase  speed  cx  is  invariant  for  the  Lindzen 
[1981]  solution,  it  varies  in  the  full  ray  solution  since  k  is  nearly 
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constant,  but  w  varies  according  to  (2),  or  equivalently  in  this 
case,  according  to  the  solution  (19).  A  region  of  positive  tidal 
shear  (Uz  >  0),  which  refracts  waves  to  a  larger  wavenumber 
\m  |  through  (1)  (since  m  <  0),  is  also  a  region  of  positive  tidal 
acceleration  (U,  >  0),  which  refracts  waves  to  higher  frequen¬ 
cies  and  larger  cx  values  through  (2).  The  two  effects  have 
opposing  influences  on  <o+  and  combine  to  reduce  the  range  of 
refraction  of  w'h  and  A,  compared  to  the  time-invariant 
Lindzen  solution.  However,  one  should  note  that  if  tidal  phase 
moved  upward,  then  U,  <  0  when  Uz  >  0,  whereupon  the  full 
ray  solution  undergoes  greater  variations  in  tu+  and  Az  than  the 
Lindzen  solution. 

As  an  extreme  example  of  the  importance  of  tidal  acceler¬ 
ations  on  gravity  wave  paths,  consider  allowing  Ar  —*  so 
that  dm/clt  — >  0  according  to  (1).  Thus  m  remains  constant 
and,  since  Kh  is  also  nearly  constant,  then  the  intrinsic  fre¬ 
quency  oT  remains  constant  through  the  dispersion  relation 
(20),  despite  the  presence  of  time-varying  background  winds. 
This  constancy  is  achieved  through  refraction  of  w  values 
through  (2),  which  in  this  case  has  the  solution  u>(f)  =  &>(/(,)  + 
k[U(t)  -  t/(f0)].  Since  w'1  =  u>  -  kU,  then  w+  remains 
constant  [e.g.,  Jones,  1969;  Zhong  et  al.,  1995]. 

The  reduced  variations  in  oj+  due  to  the  temporal  oscillation 
of  the  tide  (Figure  2c)  also  influence  the  stability  and  satura¬ 
tion  of  the  gravity  wave,  since  in  these  simulations  the  wave 
saturates  when 
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(22) 


(assuming/2  «  k2).  Since  k  remains  constant,  the  saturation 
threshold  in  (22)  varies  according  to  the  o>+  variations.  Figure 
2f  plots  the  variations  in  the  intrinsic  horizontal  phase  speed 
| cx  —  U\  —  \o)+ /k\.  We  see  again  that  waves  governed  by  the 
Lindzen  solution  are  refracted  to  smaller  intrinsic  phase 
speeds,  so  they  will  be  rendered  unstable  at  much  smaller  wave 
amplitudes  than  the  full  ray  solutions,  according  to  (22).  Con¬ 
sequently,  dissipation  is  more  likely  in  the  Lindzen  solutions. 

Figure  2g  shows  the  gravity  wave  action  density  A ,  evaluated 
in  the  absence  of  any  dissipative  processes.  We  note  that  the 
enhancements  of  A  in  the  full  ray  solution,  at  tidal  phases 
where  cg,  slows  due  to  refraction,  are  much  smaller  than  in  the 
Lindzen  solution.  This  effect  also  makes  the  waves  governed  by 
a  Lindzen  solution  much  more  susceptible  to  saturation,  since 
their  action  densities  and  associated  velocity  variances  become 
much  larger  in  response  to  refraction  than  those  governed  by 
the  full  ray  solution.  The  combined  results  in  Figures  2f  and  2g 
indicate  that  waves  governed  by  a  Lindzen  solution  are  ren¬ 
dered  unstable  more  readily  and  more  suddenly  by  tidal  re¬ 
fraction  than  similar  waves  governed  by  the  full  ray  solutions. 
Clearly,  the  time  variations  of  the  tidal  oscillation  have  a  sta¬ 
bilizing  influence  on  the  GW-tidal  interaction,  so  Lindzen  so¬ 
lutions  overestimate  tidally  modulated  gravity  wave  breaking. 

Figure  2h  shows  the  results  of  applying  saturation  limits  to 
the  action  densities  of  each  ray  solution.  As  expected,  the 
Lindzen  solution  is  more  suddenly  and  dramatically  dissipated 
at  —74  km.  To  quantify  the  resultant  mean-flow  accelerations, 
we  first  evaluate  the  vertical  flux  of  zona!  pseudomomentum 
density, 

F  =  kcgzA .  (23) 

It  can  be  shown  that  —F  in  (23)  equals  the  vertical  component 
of  the  EP  flux,  the  divergence  of  which  determines  the  accel¬ 


eration  of  the  zonally  averaged  background  flow  [e.g.,  Andrews 
and  McIntyre,  1976].  Since  here  the  EP  flux  divergence  simpli¬ 
fies  to  ~  dF/dz,  we  have  approximated  it  in  Figure  2i  by  dif¬ 
ferencing  successive  F  values  along  the  ray  path.  An  exact 
calculation  of  -dF/dz  again  requires  a  ray  tube  calculation; 
we  shall  show  later  that  this  approximate  along-ray  calculation 
is  reasonably  accurate  in  this  case. 

The  Lindzen  solution  exhibits  a  strong  and  narrow  layer  of 
EP  flux  divergence  and  mean-flow  acceleration  (— p^\z)i)F/Sz, 
shown  in  Figure  2j)  at  —74  km  and  no  response  elsewhere. 
Conversely,  the  full  ray  solution  shows  a  smaller  impulsive 
mean-flow  acceleration  due  to  saturation  at  —77  km.  However, 
this  is  superimposed  upon  a  smoothly  varying  oscillatory  mean- 
flow  acceleration,  which  arises  in  spite  of  an  absence  of  any 
dissipation  of  A  away  from  the  narrow  layer  at  —77  km.  This 
latter  feature  arises  due  to  the  transience  in  the  problem  and 
is  addressed  in  the  following  section. 


Transient  EP  Flux  Effects 

The  EP  flux  divergences  that  arise  in  the  full  ray  solution 
away  from  regions  of  wave  saturation  in  Figure  2i  can  be 
explained  by  first  multiplying  the  wave-action  equation  (8)  by 
k.  Since  horizontal  gradients  are  zero  here  and  k  is  invariant 
with  height,  this  yields 

9  clF 

—  (kA)  =  ~  (24) 


in  the  absence  of  wave  dissipation  (r =  0).  Equation  (24)  has 
the  same  form  as  the  generalized  EP  theorems  sought  in  var¬ 
ious  three-dimensional  studies  of  wave  phenomena  [e.g., 
Haynes,  1988;  Scinocca  and  Shepherd,  1992],  Nonacceleration 
conditions  require  dF/dz  =  0,  but  (24)  shows  that  the  tran¬ 
sience  in  the  wave  pseudomomentum  density  induced  by  the 
tidal  oscillation  [d(kA)/dt  0]  yields  dF/dz  +  0  and  a 
mean-flow  forcing. 

This  mean-flow  forcing  due  to  transient  pseudomomentum 
densities  differs  from  previous  studies  of  these  effects  [e.g., 
Jones  and  Houghton,  \911;  Andrews,  1980;  Coy  and  Hitchman, 
1984;  Dunkerton  and  Fritts,  1984]  in  that  the  transience  is  not 
preimposed  but  is  periodically  induced  by  tidal  refraction.  This 
leads  to  gravity  wave  EP  flux  divergences  which  are  highly 
phase  coherent  with  the  tidal  oscillation  (as  illustrated  later; 
see  Figure  3a),  which  raises  some  important  questions  con¬ 
cerning  their  possible  effects  on  tidal  wind  profiles.  We  inves¬ 
tigate  this  further  here,  keeping  in  mind  the  likely  limitations 
of  our  noninteractive  small-amplitude  slowly  varying  ray  model 
in  fully  specifying  transient  wave  mean-flow  interactions  [Fritts 
and  Dunkerton,  1984;  Fritts  et  al.,  1996], 

We  first  express  the  total  zonal  momentum  balance  of  the 
larger-scale  flow  as 


d(U  +  TIVidf.)  f  fz  +  1  1  |  ^  + 

dt  2  pa  dx  po 


dF 

dz 


(25) 


Here  (U,  V)  is  the  zonally  averaged  wind  profile,  and  D  is  the 
mean  zonal  body  force  which  sustains  the  mean  winds.  We 
assume  for  simplicity  that  D  results  from  saturation  of  small- 
scale  gravity  waves  and  does  not  influence  the  tidal  profile.  The 
extra  EP  flux  divergence  term  in  (25)  results  from  tidally  in¬ 
duced  transience  of  a  nondissipating  gravity  wave.  We  assume 
that  the  same  wave  is  launched  with  the  same  amplitude  at 
equispaced  longitudes  around  a  latitude  circle.  In  this  case, 
dF/dz  averages  to  0  over  one  horizontal  oscillation  of  the  tide 
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(a) 


(b) 


Figure  3.  Variations  with  height  of  various  GW  and  tidal 
parameters,  as  encountered  along  the  gravity  wave  path  in 
Figures  2a  and  2b:  (a)  zonal  accelerations  due  to  the  tide 
(dU'TlDE/dl,  dotted  curve)  and  transient  GW-induced  flow 
accelerations  -pg'dF/dz  (exact  solution  (27),  thick  solid 
curve,  and  the  along-ray  approximation,  thin  curve);  (b)  tidal 
winds  along  the  ray  path:  original  tidal  velocity  (/Vide  (solid 
curve)  and  modified  velocity  wave  allV,DE  (dotted  curve);  and 
(c)  upward  flux  of  zonal  pseudomomentum  density  F  for  the 
full  ray  solution  (solid  curve)  and  for  the  Lindzen  solution 
(dashed  curve).  At  z  —  60  km,  u71  +  v 71  =  1  m2  s-2. 


and  does  not  contribute  to  D.  For  generality,  we  have  assumed 
a  tidal  velocity  wave  °UVide  and  pressure  wave  S' Vide  which 
can  differ  from  the  original  values  (/Vide  and  /'Vide  as  a 
result  of  this  transient  EP  flux  term.  Note  too  that  we  assume 
K  ¥=  0  in  (25)  and  we  have  omitted  the  advective  term 
Ud/dx^’-y ide)  since  if  is  usually  small. 

Strictly,  (25)  only  holds  for  zonally  averaged  flows:  A  full 
three-dimensional  wave  mean-flow  interaction  theory  does  not 
currently  exist  [Scinocca  and  Shepherd ,  1992],  However,  the 
applicability  of  (25)  in  three  dimensions  is  fundamental  to 
nearly  all  current  specifications  of  GW-drag  effects  on  three- 
dimensional  atmospheric  circulations  [see,  e.g.,  Holton,  1984; 
Smith,  1996], 

Factoring  out  the  zonally  averaged  terms  in  (25)  yields  the 
tidal  momentum  balance  equation 


where  for  simplicity,  we  have  neglected  any  meridional  velocity 
response  of  the  tide  to  the  transient  gravity  wave  forcing.  It 
should  be  noted  that  strictly,  another  dissipation  term  should 
be  included  on  the  right  of  (26)  to  account  for  the  constant 
amplitude  with  height  of  the  original  tidal  oscillation  (rather 
than  an  exponential  increase  °cp,7  ly2(z)  that  would  occur  in 
the  absence  of  some  preexisting  tidal  dissipation). 

Thus  (26)  implies  that  these  transient  gravity  wave  EP  flux 
divergences  affect  the  tide  and  not  the  mean  flow,  consistent 
with  previous  nontransient  saturation  studies  of  GW-tidal  in¬ 
teractions  [e.g.,  Fritts  and  Vincent,  1987].  These  EP  flux  diver¬ 
gences,  approximated  in  Figure  2i  by  along-ray  differencing, 
can  be  evaluated  exactly  in  the  case  of  constant  tidal  amplitude 
and  no  gravity  wave  saturation  by  manipulating  (24),  based  on 
our  analytical  action  solution  (21),  to  give 

dF  _  BA 

dz  dt  ’ 


which  can  be  evaluated  exactly  with  the  aid  of  our  ray-tracing 
method  for  computing  dCyJdz,  as  described  in  Appendix  A. 

To  study  this,  an  identical  ray  trace  to  that  in  Figure  2  was 
performed,  except  that  in  this  case  the  gravity  wave  had  an 
initial  variance  (V77  +  i T1  -  1  m2  s~2  at  60  km.  This  smaller 
initial  variance  results  in  no  saturation  of  the  full  ray  solution 
between  60  and  100  km,  so  that  the  EP  flux  divergences  of  this 
gravity  wave  solution  are  due  entirely  to  the  transience  induced 
by  tidal  accelerations.  Figure  3a  shows  the  zonal  accelerations 
of  the  original  tidal  oscillation,  dU'-nDJdt  (dotted  curve),  as 
computed  along  the  ray  trajectory.  Zonal  accelerations  due  to 
transient  gravity  wave  EP  flux  divergences  induced  by  this  tidal 
oscillation  are  also  plotted,  using  both  the  exact  solution  (27) 
(solid  curve)  and  an  approximation  obtained  by  along-ray  dif¬ 
ferencing  (thin  curve).  Note  that  the  differencing  calculation 
slightly  overestimates  -  po'dF/dz. 

We  see  that  dll'^^dt  and  —  pg'dF/dz  =  pg  ld(kA)/dt 
have  a  similar  vertical  wavelength  and  are  essentially  in  phase, 
although  -po'dF/dz  is  clearly  nonsinusoidal  and  peaks  at 
different  heights.  This  in-phase  relationship  tends  to  enhance 
the  total  diurnal  accelerations  of  the  flow.  To  quantify  this,  we 
first  rearrange  (26)  with  the  aid  of  (24)  to  give  us 


d 

Tt 


kA\ 

^TIDE  ~  I 
Po  / 


1 

Po  dx 


(28) 


On  invoking  (28)  at  times  r,  and  t2,  during  which  we  assume 
that  O^tide/3*  's  constant,  then  it  follows  that 


^TIDeC*)  h )  ~  zi  fi) 


k\A{x,  z,  tz)  -A(x.  z,  f,)] 

PoU) 

We  now  assume  that  C/tideC*>  0  returns  after  this  unsat¬ 
urated  wave  activity  has  propagated  through  the  tide  (A  — »  0) 
[e.g.,  Andrews,  1980],  even  though  no  smooth  transition  to  A  = 
0  occurs  in  the  present  problem.  This  yields  [e.g.,  Andrews, 
1980;  Dunkerton,  1981;  Dunkerton  and  Fritts ,  1984] 


d^TIDE  1  ^9'tide_  1 

dt  po  dx  Po  dz  ’ 


(26) 


^TIDeC^i  f)  ~  (/tIDe(-*->  z>  0  + 


kA(x,  z ,  t) 
Pn(z ) 


(30) 
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and  AStide Jdx  =  0P'TlDE/dx.  Thus  we  see  that  time-varying 
gravity  wave  pseudomomentum  densities  due  to  refraction  by 
the  tidal  wind  profile  produce  accompanying  accelerations  or 
decelerations  of  the  tidal  wind  profile  which  are  temporary  and 
reversible. 

Figure  3b  shows  fy'ViDt.  as  a  function  of  height,  as  encoun¬ 
tered  along  the  ray  trajectory  in  Figure  2a.  The  dotted  curve 
shows  TlxIDE,  evaluated  from  (30),  illustrating  that  the  tran¬ 
sient  p^ld(kA)/dt  term  is  large  at  the  tidal  peak  and  tends  to 
enhance  it.  The  similar  mean-flow  acceleration  at  the  tidal 
trough  at  ~95  km  in  Figure  3b  is  due  to  the  amplification  of  a 
smaller  pseudomomentum  density  kA  by  the  smaller  back¬ 
ground  densities  p0(z)  at  these  heights,  according  to  (30). 

To  show  these  effects  more  clearly,  Figure  4  plots  time 
variations  of  the  action  density  A,  tidal  velocities  L/'T1DE  ar|d 
°^tide>  and  pseudomomentum  flux  F  at  a  fixed  height  of  80 
km.  Evaluating  time  variations  at  a  given  height  requires  that 
we  launch  a  number  of  rays  at  successive  time  intervals.  How¬ 
ever,  the  results  in  Figure  4  were  produced  by  extrapolating  the 
results  of  a  single  ray  integration  by  noting  that  for  a  constant 
tidal  amplitude,  ray  paths  at  later  launch  times  follow  by  trans- 


(a) 


(b) 


t  (hours) 


Figure  4.  Time  variations  of  GW  and  tidal  parameters  at  z  - 
80  km  for  a  continuous  train  of  gravity  waves:  (a)  action 
densities  A;  (b)  tidal  velocity  wave  LI  Vide  (solid  curve)  and 
modified  velocity  wave  °UVide  (dotted  curve);  and  (c)  zonal 
GW  pseudomomentum  flux  densities  F. 


(a) 


Figure  5.  Simulation  similar  to  that  in  Figure  4  but  using  two 
oppositely  propagating  gravity  waves  with  phase  speeds  (cx), 
and  (cx)2.  (a)  Time  variations  of  (cx),  (dashed-dotted  curve) 
and  ( cx)2  (dashed  curve)  are  plotted  at  z  =  80  km.  Tidal 
velocity  wave  U'r  1DE  (solid  curve)  and  modified  velocity  wave 
TIVioe  (dotted  curve)  are  also  shown,  (b)  Time  variations  of  F 
for  each  wave  (F,  and  Fz)  and  their  sum  FtotilI. 


lating  the  ray  loci  of  the  earlier  solution  in  Figure  2a  along  the 
phase  motion  of  the  tide  to  later  times.  Figure  4b  again  shows, 
more  clearly  this  time,  that  the  tidal  peaks  are  enhanced  by  the 
transient  gravity  wave  forcing. 

Note  that  if  we  were  to  launch  this  same  gravity  wave  with  an 
oppositely  directed  phase  speed  (cx  =  —25  m  s_l  initially),  its 
tidally  induced  transience  will  enhance  the  troughs  of  the  tidal 
velocity  wave.  Clearly,  then,  both  waves  acting  in  concert  tend 
to  enhance  both  the  peaks  and  the  troughs  of  the  tide  and  thus 
amplify  the  local  tidal  wind  oscillation.  The  results  of  such  a 
two-wave  simulation  are  shown  in  Figure  5:  Figure  5a  shows 
U' nnE  and  °U-'nnE  at  80  km,  the  latter  resulting  from  the 
combined  transience  of  two  gravity  waves  with  initial  cx  values 
of  (O,  =  25  m  s_1  (the  original  ray  solution)  and  ( cx)2  = 
-25  m  s~‘.  The  (cx)2  wave  was  launched  at  40  km  with  an 
appropriately  scaled-down  initial  amplitude,  so  that  both  waves 
had  the  same  initial  intrinsic  phase  speed  and  equal  and  op¬ 
posite  diurnally  averaged  EP  fluxes  at  each  height.  The  time 
variations  of  ( cx ) ,  and  ( cx)2  due  to  tidal  accelerations  are  also 
plotted  in  Figure  5a.  Note  particularly  the  overall  enhance¬ 
ment  of  the  local  tidal  wind  oscillation  due  to  the  combined 
transient  variations  of  both  gravity  wave  EP  fluxes.  However, 
note  again  that  these  enhancements  are  temporary;  whether 
permanent  changes  in  tidal  amplitudes  are  possible  is  dis¬ 
cussed  later. 

We  conclude  here  by  considering  whether  there  is  any  ob¬ 
servational  support  for  the  transient  gravity  wave  dynamics  we 
have  described  here.  Figures  3c,  4c,  and  5b  plot  the  pseudo¬ 
momentum  flux  densities  F.  On  comparing  these  curves  with 
the  tidal  velocity  waves  in  Figure  3b,  4b,  and  5a,  respectively, 
we  note  clear  anticorrelations  between  F  and  the  tidal  velocity 
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f/ViDE-  Since  there  was  no  saturation  in  these  particular  ray 
integrations,  these  anticorrelations  arise  solely  from  the  tran¬ 
sience  in  the  gravity  wave  amplitudes  induced  by  the  tidal 
accelerations  (for  example,  Ftotal  =  F,  +  F2  in  Figure  5b  is 
uniformly  zero  using  nondissipating  Lindzen  solutions).  This 
anticorrelation  is  a  robust  result;  it  occurs  in  varying  degrees 
within  all  of  our  simulations.  Note  particularly  from  Figure  5b 
that  the  anticorrelations  occur  regardless  of  the  direction  of 
horizontal  phase  propagation  of  the  gravity  wave. 

The  gravity  wave  EP  flux  —  F  —  for  f2  <&'  o/  <K  N2, 

where  u'  and  w'  are  zonal  and  vertical  velocity  perturbations, 
respectively,  of  the  gravity  wave  motion  [e.g.,  Fritts  and  Vincent, 
1987].  Since  the  covariance  u'  w  ’  can  be  measured,  this  permits 
observations  of  F  within  tidal  velocity  fields.  Mesospheric  ra¬ 
dar  measurements  have  revealed  that  anticorrelations  between 
F  and  U’T It>E  are  a  characteristic  feature  of  the  data  [see.  e.g., 
Wang  and  Fritts,  1991].  To  date,  such  observations  have  been 
explained  in  terms  of  tidally  modulated  gravity  wave  saturation 
[e.g.,  Fritts  and  Vincent,  1987;  Lu  and  Fritts ,  1993;  Fritts,  1995]. 
Our  results  suggest  that  some  of  this  observed  variability  may 
also  be  due  to  gravity  wave  transience  which  is  induced  by  the 
time-varying  refraction  of  gravity  waves  by  tidal  wind  oscilla¬ 
tions.  We  note  again  that  the  mechanism  we  have  described 
here  will  operate  regardless  of  whether  the  wave  field  is  dissi¬ 
pating  or  stable  and  regardless  of  the  horizontal  propagation 
directions  of  the  waves. 


Numerical  Results  for  Height-Varying 
Tidal  Amplitudes 

We  now  study  the  interaction  of  a  gravity  wave  with  a  more 
complicated  tidal  structure.  We  adopt  a  Gaussian  vertical  pro¬ 
file  of  tidal  amplitude  in  (7)  of  the  form 

tWU)  =  Um  exp  (  ~U~Zmax>)  .  (31) 

This  time  we  consider  a  semidiurnal  tide  (n  =  2)  of  Az  = 
2t t/M  =  100  km  and  amplitude  governed  by  (31),  with  f/max  = 
20  m  s'"1,  zmax  =  80  km,  and  L  =  30  km.  This  yields  tidal 
variations  in  reasonable  agreement  with  mesospheric  data  dur¬ 
ing  summer  at  35°S  [e.g.,  Vincent  et  at.,  1988],  although  one 
should  note  that  significant  semidiurnal  tidal  activity  often 
persists  well  into  the  thermosphere  [e.g.,  McLandress  et  ai, 
1996].  These  tidal  winds  are  contoured  in  Figure  6a. 

The  transformation  sequence  (x,  z)  -»•  ( x ",  z")  now  no 
longer  eliminates  time  variations  in  the  (x",  z”)  frame,  since 
the  amplitude  profile  (31)  moves  upward  with  time  when 
viewed  from  this  frame.  Consequently,  we  must  now  resort 
solely  to  numerical  ray  integrations  to  obtain  gravity  wave 
solutions. 

We  start  the  rays  in  these  simulations  from  an  artificial  lower 
height  of  z  =  —50  km,  where  U(z,  t)  and  its  height-time 
gradients  are  all  extremely  small.  We  launch  rays  of  identical 
characteristics,  so  that  rays  emanating  from  well  below  the  tide 
are  initially  monochromatic  and  plane-parallel  ( dm/dz  =  0). 
In  this  way  the  ray  integrations  from  z  =  —50  km  to  the 
surface  are  a  numerical  technique  for  establishing  accurate 
initial  conditions  for  the  wave  at  the  ground  under  these  as¬ 
sumptions. 

Figure  6  shows  the  results  of  a  numerical  ray  integration, 
with  variations  of  parameters  along  the  ray  plotted  in  the 
height  range  0-100  km.  The  wave  had  identical  original  pa¬ 


rameters  to  the  one  traced  in  Figure  2,  except  that  it  was 
launched  at  t  =  0  and  z  =  —  50  km  with  an  initial  u’1  +  v'~ 
that  grows  to  7  m2  s-2  at  z  =  60  km  in  the  absence  of  any 
refraction.  This  choice  ensures  that  the  wave-action  densities 
here  are  of  similar  magnitude  to  those  in  Figure  2. 

Below  ~40  km,  where  U(z,  t)  is  small,  we  note  that  the 
Lindzen  and  full  ray  solutions  agree.  Above  this  height,  where 
the  time-varying  tidal  amplitudes  become  significant,  we  note 
again  that  the  two  solutions  diverge  significantly  from  each 
other. 

To  study  the  behavior  of  the  full  ray  solutions  in  greater 
depth,  Figure  7a  shows  time-height  ray  loci  for  a  sequence  of 
successive  ray  integrations  for  this  wave  that  commenced  at 
z  =  -50  km  at  regular  15-min  intervals.  Again,  the  tidal 
velocity  oscillation  through  which  the  waves  propagate  is  con¬ 
toured  on  the  plot.  This  time  the  waves  are  traced  to  a  height 
of  150  km,  where  the  tidal  oscillation  is  almost  completely 
absent.  This  is  done  purely  to  provide  a  “far  field”  from  the 
tide  which  allows  us  to  gauge  whether  interaction  with  the  tide 
permanently  affects  the  emergent  wave  activity.  Consequently, 
no  attempt  has  been  made  to  include  ionospheric  processes 
into  the  background  atmosphere  at  heights  >  100  km. 

The  ray  loci  in  Figures  7a  and  7b  reveal  that  permanent 
changes  in  the  ray  paths  result  from  the  interaction  with  the 
tide,  so  that  at  —150  km  there  is  a  periodic  (semidiurnal) 
focusing  and  defocusing  of  ray  paths  that  were  all  initially 
parallel.  While  the  ray  focusing  events  in  Figure  7a  look  like 
possible  caustics,  they  are  not,  as  evident  in  Figure  7b,  while 
Figure  7f  shows  that  WKB  requirements  are  satisfied.  The 
focusing  of  rays  implies  a  compression  of  wave  packets  and  a 
consequent  buildup  of  wave  amplitude,  whereas  reduced  am¬ 
plitudes  should  occur  during  the  defocusing  phase  as  packet 
widths  broaden.  The  action  solutions  in  Figure  7c  bear  this  out. 
Note  again  that  these  time-varying  action  densities  result  solely 
from  time-varying  tidal  accelerations  acting  on  the  ray  paths. 
The  horizontal  line  in  Figure  7c  shows  the  time-mean  action 
density  of  all  the  rays  at  —150  km.  which  also  equals  the 
uniform  A  values  that  occur  here  when  Lindzen  solutions  arc 
used  to  evaluate  the  ray  paths.  Figure  7e  shows  a  similar 
collection  of  curves  for  o>+. 

The  time-varying  ray  characteristics  at  150  km  occur  because 
the  packet  width  of  the  tide  is  of  the  order  of  the  vertical  tidal 
wavelength  (LI A,  —  I).  When  L/Az  »  1,  the  time  varia¬ 
tions  of  the  emergent  wave  activity  are  substantially  reduced 
because,  after  propagating  through  many  tidal  cycles  during 
which  tidal  amplitudes  change  only  slowly,  the  time-varying 
refraction  terms  tend  to  average  out  to  a  near-zero  net  effect 
on  the  gravity  wave. 

The  time-varying  A  values  at  150  km  imply  transient  mean- 
flow  forcing  of  a  form  proportional  to  HA  lilt,  which  is  shown  in 
Figure  7d.  The  buildup  of  action  densities  during  focusing 
phases  also  suggests  that  fairly  intense  wave  breaking  will  occur 
at  times  of  ray  focusing.  Since  the  focusing  is  periodic,  such 
wave  breaking  may  force  a  wave  motion,  particularly  in  the 
presence  of  a  second  oppositely  propagating  gravity  wave  as  in 
Figure  5.  On  inspecting  the  upward  sloping  time-height  focus¬ 
ing  regimes  in  Figure  7a,  wave  breaking  within  these  regions 
would  tend  to  produce  a  semidiurnal  wave  with  upward  phase 
propagation  and  a  very  long  vertical  wavelength.  If  this  pro.- 
duced  an  internal  gravity  wave  mode,  this  wave  would  propa¬ 
gate  to  lower  heights.  However,  since  the  forcing  is  also  highly 
nonsinusoidal  (Figure  7d),  the  response  to  this  forcing  may  be 
more  complicated. 


ECKERMANN  AND  MARKS:  RAY  MODEL  OF  GRAVITY  WAVE-TIDAL  INTERACTIONS 


21,205 


(a) 


(c) 


(g) 


1 


too 
80 
60 
40 
20 
0 

-0.010  0.000  0.010  0.020  0.030  0.040 


-dF/dz  (kg  m'2  s’1  day ‘) 


(b) 


(f) 


-pQ  dF/dz  (m  s’1  day1) 


Figure  6.  As  for  Figure  2,  but  for  a  semidiurnal  tide  with  Az  —  2 irlM  =  100  km  and  a  Gaussian  amplitude 
profile  governed  by  (31 ),  with  L/nlIlx  =  20  ms1,  zm.,x  =  80  km,  and  L  =  30  km.  Results  are  plotted  from  the 
ground  to  z  =  100  km.  (a)  Ray  loci  through  semidiurnal  tidal  winds,  (b)  spatial  trajectory,  (c)  intrinsic 
frequency,  (d)  vertical  wavelength,  (e)  ground-based  zonal  phase  speed,  (f)  intrinsic  zonal  phase  speed,  (g) 
unsaturated  action  density,  (h)  saturation-scaled  action  density,  (i)  Eliassen-Palm  flux  divergence,  and  (j) 
mean-flow  acceleration. 


Discussion 

Lindzen  Parameterizations  of  GW-Tidal  Interactions 

The  ray  simulations  have  illustrated  that  Lindzen  solutions 
often  differ  markedly  from  the  full  ray  solutions  within  a  time- 
varying  tidal  background.  This  raises  some  doubts  about  the 
accuracy  of  the  Lindzen  [1981]  parameterization  (and  subse¬ 
quent  variants  of  it)  in  modeling  GW-tidal  interactions. 

The  Lindzen  parameterization  remains  popular  among 
modelers  because  of  its  simplicity  and  its  generally  good  per¬ 
formance  in  simulating  mean  circulation  profiles  within  arbi¬ 
trary  environments.  This  contrasts  with  the  simulations  con¬ 
ducted  here,  which  needed  a  fairly  complicated  numerical 
model  to  accurately  solve  a  highly  idealized  problem.  Thus  we 


consider  here  whether  current  implementations  of  the  Lindzen 
parameterization  can  be  modified,  rather  than  overhauled  or 
replaced,  to  include  some  of  the  time-varying  effects  we  have 
reported  here,  while  retaining  some  of  the  attractive  simplicity 
of  the  current  parameterization. 

Of  the  many  changes  to  gravity  wave  solutions  that  arise  on 
including  background  time  variations  there  are  two  major  ones 
which  have  important  consequences  for  the  Lindzen  [1981] 
parameterization.  The  first  is  the  different  solution  of  the 
wave-action  equation  (8).  The  Lindzen  GW-drag  formulae  are 
derived  from  the  solution  (12)  [e.g.,  Schoeberl,  1985],  which 
follows  from  (8)  on  neglecting  ,v,  y,  and  t  variations  in  the 
background  atmosphere.  Retention  of  time-varying  terms  re- 
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Figure  7.  Mulliray  simulation  through  the  same  tidal  profile  as  in  Figure  6.  (a)  Zonal  tidal  winds  are 
contoured  from  the  ground  to  150  km;  thin  curves  show  various  ray  loci  through  these  semidiurnal  tidal  winds, 
(h)  horizontal  trajectories;  x  —  t  loci  of  the  various  rays  traced  through  the  tide,  (c)  action  densities,  A ,  (d) 
rate  of  change  of  A ,  dA/dt,  (e)  intrinsic  frequencies,  to+,  and  (f)  WKB  parameters  8  and  y.  Figures  7c-7f  show 
time  variations  of  the  ray  parameters  at  z  —  150  km. 


quires  a  numerical  solution  of  (9),  as  described  in  Appendix  A. 
Although  analytical  solutions  to  (9)  can  be  derived  for  certain 
problems,  such  as  (16)  in  the  case  of  a  constant  tidal  amplitude, 
even  these  solutions  have  possible  ramifications  for  the  gravity 
wave  driving  of  the  larger-scale  flow  which  are  difficult  to 
delineate  fully  at  present.  As  any  general  implementation  of 
(9)  would  require  a  major  overhaul  of  the  Lindzen  parameter¬ 
ization,  we  investigate  here  whether  improvements  can  result 
without  changing  the  current  GW-drag  formulae  based  on 
(12). 

A  second  major  effect  of  time  variations  is  that  ground-based 
horizontal  phase  speeds  c.x  vary  in  response  to  background 
accelerations,  as  shown  in  Figures  2e  and  6e  [see  also  Jones  and 
Houghton,  1971;  Frills  and  Dunkerlon,  1984;  Walterscheid  and 
Schubert ,  1990;  Zhong  et  al.,  1995].  By  contrast,  the  Lindzen 


theory  assumes  no  background  time  variations,  socr  is  taken  to 
be  constant  with  height.  Variations  in  cx  are  important  because 
the  parameterization  formulae  of  Lindzen  [1981]  are  especially 
sensitive  to  phase  speed  values;  for  example,  the  turbulent 
diffusivities  that  result  from  wave  breaking  scale  as  |cr  —  U |4. 

Representative  variations  in  cx  can  be  incorporated  fairly 
easily  into  the  Lindzen  parameterization.  One  simple  possible 
implementation  is  described  and  discussed  in  Appendix  B.  To 
see  whether  including  transient  cx  variations  alone  leads  to  a 
better  comparison  with  a  full  ray  solution.  Figure  8  reproduces 
the  simulation  of  Figure  2  but  includes  an  extra  solution  (the 
thick  dotted  curve)  in  which  the  full  trajectory  solution  is  used 
(i.e.,  cx  can  vary),  but  wave  amplitudes  are  governed  by  the 
time-invariant  solution  (12).  Clearly,  the  inclusion  of  variable 
cx  values  improves  the  simulation.  Thus  an  adapted  Lindzen 
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Figure  8.  As  for  Figure  2,  but  including  a  third  solution  in  which  the  full  ray  solution  is  used  for  the  ray  path 
(so  that  cx  variations  are  included).  However,  the  Lindzcn-like  solution  (12)  is  used  to  govern  the  wave 
amplitude.  This  “adapted  Lindzen”  solution  is  plotted  with  a  thick  dotted  curve,  in  addition  to  the  original  full 
ray  solution  (solid  curve)  and  Lindzen  solution  (dashed  curve),  (a)  Ray  loci  through  diurnal  tidal  winds,  (b) 
spatial  trajectory,  (c)  intrinsic  frequency,  (d)  vertical  wavelength,  (e)  ground-based  zonal  phase  speed,  (f) 
intrinsic  zonal  phase  speed,  (g)  unsaturated  action  density,  (h)  saturation-sealed  action  density,  (i)  Eliasscn- 
Palm  flux  divergence,  and  (j)  mean-flow  acceleration. 


parameterization,  in  which  the  original  amplitude  formulae  are 
retained  but  cx  values  may  vary,  should  result  in  more  accurate 
simulations  of  gravity-wave  interactions  with  tides. 

Transient  GW-Induced  Tidal  Accelerations 

The  ray-tracing  experiments  simulated  time-varying  EP  flux 
divergences  for  a  gravity  wave  as  it  propagated  through  a  tide. 
It  was  shown  that  this  produced  zonal  accelerations  which  were 
highly  phase  coherent  with  the  tide,  and  it  was  argued  that  they 
temporarily  amplified  the  tide. 

For  permanent  tidal  amplification  or  attenuation  to  occur 
the  total  energy  density  of  the  gravity  wave  packet  must  be 
changed  by  the  interaction.  If  A  behaves  conservatively  (r- 1  = 


0),  then  this  requires  a  net  change  in  «+  [e.g.,  Broutman  and 
Young,  1986].  Broutman  and  Young  [1986]  and  Broutman  and 
Grimshaw  [1988]  found  large  permanent  increases  in  w+  for 
oceanic  gravity  waves  propagating  through  an  inertial  oscilla¬ 
tion.  but  these  required  c„2  =  Cz  somewhere  along  the  ray 
path,  which  cannot  occur  in  these  simulations.  For  permanent 
changes  in  tu+  to  occur  here,  transient  GW-induced  modifica¬ 
tions  to  the  tide  must  feed  back  to  “self-accelcrate”  later  por¬ 
tions  of  the  wave  packet  that  propagate  into  this  region.  Fig¬ 
ures  3a,  3b,  and  4b  show  that  asymmetric  °lf  j  [DL;,  SPVJnoyaz, 
and  'rurrJdl  profiles  arise  from  transient  forcing  by  a  single 
wave.  On  substituting  these  terms  into  (1)  and  (2)  we  see  that 
subsequent  packets  which  propagate  through  the  modified 
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wind  profile  ^U-Vide  can  undergo  small  net  changes  in  m  and  to 
and  hence  to+  on  propagating  through  one  full  tidal  cycle.  If 
to h  increases  (decreases),  wave  packets  which  propagate 
through  gain  (lose)  energy  during  the  interaction  and 

thus  attenuate  (amplify)  the  tide  [e.g.,  Broutman  and  Grirn- 
sliaw,  1988], 

These  arguments,  however,  ignore  saturation  and  a  host  of 
other  possible  nonlinear  influences.  To  assess  these  effects,  a 
more  complete  numerical  model  is  required  which  incorpo¬ 
rates  finite-amplitude  wave  mean-flow  interactions.  The  next 
logical  extension  from  here  would  probably  be  a  “quasi-linear” 
numerical  model  of  gravity  wave  self-acceleration  and  satura¬ 
tion  processes  within  tidal  profiles,  following  similar  models  of 
these  processes  within  mean  flows  [e.g.,  Jones  and  Houghton, 
1971,  1972;  Dunkerton  and  Fritts,  1984;  Mobhs,  1985;  Tanaka, 
1986;  Tanaka  and  Yoshizawa,  1987].  However,  comprehensive 
nonlinear  models  of  transient  gravity  wave  mean-flow  interac¬ 
tions  have  suggested  that  certain  aspects  of  the  evolution  are 
intrinsically  nonlinear  and  thus  beyond  the  scope  of  any  slowly 
varying  small-amplitude  quasi-linear  model  [Fritts  and  Dunker¬ 
ton,  1984;  Walterscheid  and  Schubert,  1990;  Fritts  et  al.,  1996], 
Thus  nonlinear  modeling  approaches  may  ultimately  be 
needed  to  specify  fully  the  role  of  gravity  wave  transience  in 
GW-tidal  interactions.  Since  tides  are  just  large-scale  gravity 
waves,  this  work  may  also  aid  modeling  of  GW-OW  interac¬ 
tions  involving  a  big  wave  and  short  wave,  which  are  believed 
to  be  important  in  the  formation  of  atmospheric  GW  spectra 
[Hines,  1991], 

Of  further  interest  were  the  numerical  simulations  in  Figure 
7  which  showed  that  the  interaction  of  a  train  of  initially  steady 
monochromatic  plane-parallel  gravity  waves  with  an  isolated 
tidal  structure  could  induce  permanent  d(kA)/dt  variations  in 
the  emergent  wave  field.  It  w7as  noted  that  at  higher  altitudes, 
wave  dissipation  at  “focusing  phases”  could  force  a  disturbance 
with  a  similar  tidal  period  but  that  it  would  probably  propagate 
downward.  These  findings  raise  some  interesting  questions 
concerning  diurnally  varying  gravity  wave  activity  and  its  pos¬ 
sible  influence  on  tides  in  the  middle  and  upper  atmosphere. 

It  has  been  argued  that  diurnal  variability  in  tropical  con¬ 
vective  zones  provides  a  strong  thermal  forcing  of  tidal  modes 
[e.g.,  Lindzen,  1978;  Hamilton,  1981;  Williams  and  Avery,  1996], 
Bergman  and  Salby  [1994]  found  that  direct  convective  forcing 
of  planetary-scale  tides  was  limited;  it  tended  to  force  diurnal 
oscillations  over  a  broad  range  of  spatial  scales.  They  also 
inferred  strong  gravity  wave  EP  flux  production  by  convection. 
Observations  have  also  indicated  that  diurnally  varying  con¬ 
vection  produces  a  diurnal  variation  in  the  intensity  of  convec- 
tively  generated  high-frequency  gravity  waves  entering  the 
middle  atmosphere  [e.g.,  Nastrom  and  Eaton ,  1995;  Sato  el  al., 
1995],  This  implies  a  diurnally  varying  d(kA)/dt  for  gravity 
waves  produced  by  convective  sources  in  the  lower  atmo¬ 
sphere.  Depending  on  its  phase  relationship  with  existing  tidal 
accelerations  i>U’Tlc>E/dt,  such  diurnal  d(kA)/dt  values  will 
either  amplify  or  attenuate  tides  in  the  middle  atmosphere, 
according  to  (26)  and  (28),  and  might  provide  an  indirect  route 
for  convective  tidal  forcing. 


Lindzen  parameterization  of  these  dynamics,  substantially  al¬ 
tered  predictions  about  both  the  locations  and  amounts  of 
gravity  wave  breaking  within  the  tide  result.  We  demonstrated 
that  correcting  the  ray  path  (thus  including  cx  variations)  with¬ 
out  correcting  the  wave  amplitude  brings  the  calculations 
closer  to  the  full  ray  solution.  Consequently,  we  proposed  an 
adapted  Lindzen  parameterization,  which  incorporated  cx  vari¬ 
ations  but  retained  its  original  nontransient  wave-amplitude 
formulae,  for  use  in  improved  modeling  of  GW-tidal  interac¬ 
tions  and  perhaps  other  GW  interactions  with  a  time-evolving 
background. 

Even  in  the  absence  of  wave  breaking,  our  amplitude  calcu¬ 
lations  (based  on  wave-action  conservation)  revealed  that 
time-varying  gravity  wave  EP  flux  divergences  resulted  from 
tidal  refraction  and,  consequently,  were  highly  phase  coherent 
with  the  tidal  oscillation.  Elementary  analysis  showed  that  a 
temporary  increase  in  the  amplitude  of  the  tidal  velocity  oscil¬ 
lation  resulted  and  that  permanent  energy  transfers  between 
the  tide  and  gravity  wave  might  occur  in  a  more  complete 
interactive  simulation  of  these  processes.  As  some  initial  ob¬ 
servational  support  for  the  transient  GW  amplitudes  that  we 
simulated,  we  demonstrated  characteristic  anticorrelations  be¬ 
tween  tidal  velocity  oscillations  and  gravity  wave  pseudomo¬ 
mentum  flux  densities  in  agreement  with  observations.  In 
short,  our  simulations  have  stressed  the  importance  of  tidal 
time  oscillations  in  obtaining  an  accurate  description  of  GW- 
tidal  interactions. 


Appendix  A:  Numerical  Ray-Tracing  Solution 
of  the  Wave-Action  Equation 


In  the  first  version  of  GROGRAT  [Marks  and  Eckermann, 
1995],  possible  time  variations  in  the  model  atmosphere  were 
not  included.  Thus  dc„z/dt  =  0,  where  cgz  —  dzldt  is  the 
ground-based  vertical  group  velocity  of  the  wave,  which  then 
enabled  the  wave-action  continuity  equation  (8)  to  be  reex¬ 
pressed  as  a  pseudo-ray-tracing  equation  for  F  —  c//zA .  Since 
time  variations  in  the  background  atmosphere  are  included 
here,  then  3cgz/dt  +  0  and  the  current  implementation  of  the 
wave-action  equation  is  no  longer  valid.  Here  we  present  a 
numerical  solution  of  the  wave-action  equation  which  is  valid 
in  atmospheres  which  vary  in  both  height  and  time. 

Our  implementation  is  similar  to  that  of  Broutman  [1984]. 
On  ignoring  background  horizontal  gradients  and  dissipative 
effects  the  wave-action  equation  can  be  expressed  in  the  pseu¬ 
do-ray-tracing  form  (9),  i.e., 


dA 


Dc„ 


dt  ^  dz 


(Al) 


The  form  is  “pseudo-ray  tracing”  because  in  general  3cgJdz 
cannot  be  evaluated  using  a  single  ray  but  instead  requires 
simultaneous  tracing  of  a  “tube”  of  closely  spaced  rays  of 
similar  initial  characteristics. 

However,  we  can  avoid  ray  tube  calculations  here  on  manip¬ 
ulating  (Al).  First,  we  express  the  dispersion  relation  (11)  in 
the  form 


Summary  and  Conclusions 

Ray  simulations  of  gravity  wave  propagation  through  an 
idealized  tidal  background  have  shown  that  tidal  accelerations 
modify  not  only  ray  paths,  as  was  shown  by  Zhong  et  al.  [  1 995], 
but  also  gravity  wave  amplitudes.  When  compared  to  a 


a >(z,  t )  =  Y(k,  z,  t),  (A2) 

where  k  =  (k,  l,  m),  (k,  z,  /)  is  the  so-called  augmented 
space,  and  (z,  t)  is  the  propagation  space  [Hayes,  1970].  The 
ray  equations  follow  from  (A2)  and  yield  [Hayes,  1970; 
Whitham,  1974] 
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cyz=  dY/dm.  (A3) 

Since  m  =  m(z,  t),  then  (A2)  and  (A3)  can  be  used  to 
reexpress  (A1)  as 

dA  I  r)2Y  dm  d2Y 

=  —  A  I - 1 - 

dt  \  dm 2  dz  dz  dm 


Since  the  mean  vertical  wind  velocity  W  —  0  here,  then  using 
the  relations  of  Marks  and  Eckermann  [  1995], 


d2Y 

dm2 


-ik'  +  r-w1-/2) 

u/A2 


(A5) 


launching  monochromatic  waves,  plane-parallel  ray  paths  re¬ 
sult  (see  Figure  7a)  for  which  dm/dz  =  0  initially. 

The  dm/dz  values  are  also  used  to  check  that  the  WKB 
approximation  remains  valid  along  the  ray  path,  using  (3). 
They  can  also  be  used  to  evaluate  (5)  through  the  relation 

0o>+  dm  1 

—  =  -kcJJz -cg!  —  + 


■  [(A2), (A'2  +  l2)  -  (<x2),(W+!  ~/2)]  - 
•  [(A2)2(*2  +  l2)  -  («WJ  -/*)].  (A9) 


d2Y  -(k2  +  l2)m 
dz  dm 


<u+A: 


3  (a2),  + 


(A2); 


(A2  -f2) 


2A 


(k2  +  P)(N\+fHa2); 


(A6) 


where  subscripted  background  atmospheric  variables  denote 
local  partial  derivatives  of  that  parameter  (i.e.,  derivatives  in 
the  propagation  space). 

We  are  now  left  with  the  problem  of  computing  dm/dz  in 
(A4),  which  generally  requires  the  ray  tube  calculation.  An 
elegant  way  of  calculating  it  within  the  backgrounds  considered 
here  is  to  use  the  “derived  ray  equations”  of  Hayes  [1970], 
following  Broutman  [1984],  On  taking  d/dz  to  both  sides  of  the 
ray  equation  dm/dt  —  -dY/dz  one  obtains  a  ray-tracing 
equation  for  dm/dz ,  which,  on  substituting  the  expression  for 
dCgJ dz  from  (A4),  takes  the  form 


d  I  dm 
dt  \  dz 


d2Y  i  dm\  a2Y 


\  a2Y  /am\2a2Yl 
L  dz 2  dz  )  dz  dm  +  \  dz  )  dm2\  ’ 


+  2 


(A7) 
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In  the  simulations  conducted  here,  background  gradients  in  V, 
A2,  and  a2  are  ignored,  whereupon  (A6)  is  identically  zero  and 
the  right-hand  side  of  (A8)  simplifies  to  kUIz. 

Thus  dm/dz  can  be  evaluated  along  a  ray  by  integrating 
(A7)  in  conjunction  with  the  other  ray  equations.  However, 
one  still  requires  its  initial  value  in  order  to  do  this.  Approxi¬ 
mations  such  as  dm/dz  «*  cgz  dm/dt  are  inadequate,  as  the 
wave-action  density  A,  evaluated  along  the  ray  using  (A4), 
proves  acutely  sensitive  to  even  slight  initial  errors  in  dm/dz. 
For  the  simulations  outlined  in  this  paper  the  following  ap¬ 
proaches  were  adopted.  For  tides  of  constant  amplitude,  grav¬ 
ity  waves  were  launched  either  at  peaks  or  troughs  of  the  tide 
for  which  dm/dz  =  0  is  an  accurate  initial  condition  when 
temperature  and  density  gradients  arc  ignored.  For  modulated 
tidal  amplitudes  we  launch  gravity  waves  of  identical  initial 
parameters  at  regularly  spaced  time  intervals  from  a  lower 
height  where  the  tidal  amplitude  was  very  close  to  zero.  On 


The  above  equation  follows  from  first  expanding  da>+ /dt  using 
the  definition  o)+  =  at  -  kU  and  then  simplifying  using  the 
ray-tracing  equations  for  o >  and  m,  the  phase  consistency  re¬ 
lation  dm/dt  =  —  dm/dz,  and  discarding  horizontal-gradient 
terms  (which  also  imply  a  constant  k).  A  relation  similar  to 
(A9)  was  derived  by  Broutman  [1984]. 

Appendix  B:  Time-Varying  Phase  Speeds  Within 
a  Lindzen  Parameterization 

Here  we  briefly  describe  a  method  whereby  time-varying  cx 
values  can  be  incorporated  into  GW-drag  parameterizations 
based  on  the  Lindzen  [1981]  formulation.  These  changes  have 
not  been  tested  in  a  large-scale  model,  so  there  may  be  better 
ways  of  coding  these  effects  operationally. 

We  envisage  an  array  of  ground-based  horizontal  phase 
speeds  cijn,  where  i  and  j  are  array  indices  at  model  points 
along  thex  andy  directions,  respectively,  and  n  is  the  “nth”  ray 
at  this  horizontal  location  (;,  j).  In  common  with  many  im¬ 
plementations  of  the  Lindzen  parameterization  we  ignore 
oblique  wave  propagation  paths  and  consider  only  the  vertical 
wave  propagation  above  the  starting  horizontal  location 
The  current  height  of  the  various  waves  is  stored  in  an  array 

Z‘,j.n- 

Tn  the  Lindzen  parameterization,  cx  values  are  constant, 
whereas  here  we  allow  them  to  change  in  response  to  time- 
varying  background  conditions.  If  we  assume  that  all  of  these 
waves  have  the  same  (constant)  zonally  aligned  horizontal 
wavenumber  k,  then  a>*j  „  =  k(Cj  —  U)  and  mi  J  n  is  given 
by  the  dispersion  relation,  where  the  background  zonal  wind 
speed  U  is  evaluated  at  the  wave’s  current  location  (as  given  by 
i,  /,  and  zi  Jt„)  and  at  the  current  time  t. 

Before  going  on  to  the  next  time  step  of  a  model  run  we 
evaluate  dm/dt  for  each  mi  ;  „  value  using  (1),  again  at  the 
position  (/,  j,  z, ,  „)  that  the  wave  now  occupies.  Since  most 
implementations  of  the  Lindzen  parameterization  are  hydro¬ 
static  and  use  zonally  aligned  waves,  only  the  first  and  third 
terms  on  the  right  of  (1)  are  needed.  The  third  term  could 
easily  be  omitted  too,  since  its  effects  are  generally  small.  If 
retained,  the  third  term  should  be  simplified  to  reflect  a  hy¬ 
drostatic  dispersion  relation. 

If  the  model  time  step  is  At,  then  mj  j  tl  is  incremented  by  an 
amount  At(dmj  j  nldt).  Next  we  evaluate  <x>Y  „  for  each  ray 
using  the  (hydrostatic)  dispersion  relation.  We  can  then  update 
the  phase  speed  values  as  c,  ;>„  =  « tj.Jk  +  U.  Finally,  we 
use  these  new  parameters  to  propagate  the  wave  packet  to  its 
new  height  by  incrementing  zi  j  n  by  At(c gz)lln. 

These  updated  parameters  can  then  be  used  to  quantify  the 
GW-drag  on  the  atmosphere  at  the  next  time  step  of  the  model 
using  the  formulae  of  the  Lindzen  parameterization  (some 
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vertical  interpolation  may  be  required  to  regrid  the  results  at 
the  irregularly  spaced  heights  zijn  onto  the  regular  height  grid 
of  the  model).  If  time  and  height  variations  of  the  background 
atmosphere  are  present  and  horizontal  variations  are  small, 
then  ci  J  n  will  now  vary  implicitly  according  to  (2).  One  could 
also  integrate  (2)  explicitly  for  each  wave  as  a  check  of  the 
numerics. 

In  simulations  of  GW-tidal  interactions  these  modifications 
should  be  fairly  robust,  as  illustrated  in  Figure  8.  The  exact 
response  of  a  complex  circulation  model  to  this  modified  pa¬ 
rameterization  is  harder  to  gauge.  For  example,  one  clear  side 
effect  is  that  gravity  wave  drag  itself  drives  SU/dt  variations, 
which  will  then  feed  back  via  (2)  to  accelerate  the  original  c;  j  n 
values  of  the  gravity  waves.  However,  such  “self-acceleration” 
effects  have  been  observed  in  model  simulations  and  may  be  a 
real  feature  of  the  gravity-wave  driving.  WKB  methods  such  as 
these  are  incapable  of  fully  simulating  all  aspects  of  these 
interactions  [Fritts  and  Dunkerton,  1984;  Mobbs,  1987;  Waltcr- 
scheid  and  Schubert ,  1990;  Fritts  et  al. ,  1996],  but  could  they 
simulate  the  first-order  effects? 

To  investigate  this,  we  look  to  the  results  of  numerical  mod¬ 
els.  Perhaps  the  most  detailed  study  to  date  of  gravity  wave 
dissipation  and  self-acceleration  within  evolving  mean  flows 
was  reported  by  Walterscheid  and  Schubert  [1990],  based  on 
experiments  with  a  nonlinear  two-dimensional  numerical 
model.  They  found  that  dissipating  gravity  waves  which  accel¬ 
erated  the  background  flow  tended  to  self-accelerate  in  a  way 
that  kept  their  intrinsic  phase  speeds  approximately  constant 
with  height,  a  process  they  were  unable  to  explain. 

Ray  theory  seems  to  provide  a  simple  explanation  of  this 
feature.  In  the  model  runs  of  Walterscheid  and  Schubert  [1990], 
r W/dz  was  generally  small,  but  dU/dt  was  large  and  positive 
due  to  strong  mean-flow  accelerations  induced  by  breaking 
gravity  wave  packets  [see  Walterscheid  and  Schubert,  1990,  Fig¬ 
ure  4],  As  successive  wave  packets  propagate  into  these  regions 
of  accelerating  wind,  if  spatial  wind  gradients  are  small  but 
dU/dt  is  large  and  is  accelerating  U  toward  the  ground-based 
phase  speed  of  the  wave  (as  is  typical  of  gravity-wave-induced 
mean-flow  driving),  then  m  remains  approximately  constant 
from  (1)  and  <o  refracts  according  to  (2)  to  keep  a>~  approxi¬ 
mately  constant  along  the  ray  [e.g.,  Jones,  1969].  Since  k  also 
remains  constant  if  horizontal  gradients  of  the  background 
wind  are  small,  then  constant  ru+  implies  an  approximately 
constant  intrinsic  phase  speed  \cx  —  U\,  as  found  by  Walters¬ 
cheid  and  Schubert  [1990], 

From  this  it  seems  possible  that  a  modified  Lindzen  param¬ 
eterization  which  permits  self-accelerated  cx  values  might  even 
improve  its  simulation  of  mean-flow  driving  processes. 
Whether  in  practice  self-accelerated  cx  values  improve  or  de¬ 
grade  the  stability  and/or  performance  of  the  parameterization 
in  a  larger-scale  model  remains  to  be  seen.  Of  the  existing 
GW-drag  specifications  the  models  of  Tanaka  [1986]  and 
Tanaka  and  Yoshizawa  [1987]  have  most  in  common  with  the 
modified  Lindzen  parameterization  described  here  and  pro¬ 
duced  mean  circulations  which  were  stable  and  basically  in  line 
with  observations. 
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